1 



Reconnection of Unstable/Stable Manifolds of the Harper Map 

Asymptotics-Beyond- All-Orders Approach 

Shigeru Ajisaka*) and Shuichi Tasaki 

Advanced Institute for Complex Systems and Department of Applied Physics, 
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The Harper map is one of the simplest chaotic systems exhibiting reconnection of invari- 
ant manifolds. The method of asymptotics beyond all orders (ABAO) is used to construct 
unstable/stable manifolds of the Harper map. By enlarging the neighborhood of a singu- 
larity, the perturbative solution of the unstable manifold is expressed as a Borel summable 
asymptotic expansion in a sector including t = —oo and is analytically continued to the other 
sectors, where the solution acquires new terms describing heteroclinic tangles. When the pa- 
rameter changes to the reconnection threshold, the unstable/stable manifolds are shown to 
acquire new oscillatory portion corresponding to the heteroclinic tangle after the reconnec- 
tion. 

§1. Introduction 

Bifurcation involving chaotic motion is one of tlie origins of diversity in com- 
plex systems and reconnection among unstable/stable manifolds is such an example. 
However, general features of the interplay between bifurcation and chaos are not well 
understood and a case study on a simple system would provide useful information 
for the generic behavior. Moreover, since the trajectory near the bifurcation point 
is considered to be very complicated, an analytical approach would be more helpful. 
Thus, as one of the simplest systems exhibiting reconnection,^'!^ we analytically 
study the nearly integrable Harper map. 

The Harper map depends on a real parameter k and is defined on the torus 
{{v,u) G [-7r,7r]2}: 

v{t + (t) — v{t) = — iTsinn(t) 
u{t -\- a) — u{t) = ka smv{t + a) 

where a{> 0) is the time step and it plays a role of the small parameter. Since the 
case of /c < is conjugate to that of /c > cP and since the solution for k > \ can be 
obtained from that for A; < 1 by a simple symmetry argument (cf. Appendix^, it 
is sufficient to consider the case of < < 1. In the continuous time limit, where 
(T — > 0, the map reduces to an integrable set of differential equations: 

V {t) = — smu{t) 

u'{t) = ksmv{t) . (1-2) 

Throughout this paper, the prime is used to indicate the differentiation with respect 
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to time t such as u'{t) = . Eq. (|1-2|) admits topologically different separatrices 
depending on the parameter k (cf. Fig. ^ and the separatrix changes its shape when 
— > 1. As shown in the figure, when k increases, the middle point of the upper 
separatrix moves upward and, when A; = 1, it reaches the fixed point (0,7r). This 
view would be changed for nonvanishing a due to the heteroclinic tangles of stable 
and unstable manifolds. Indeed, in case of < /c < 1, the unstable manifold of 
the fixed point (vr, 0) is expected to show an oscillatory behavior near (— vr, 0), while 
in case of k = 1, it would oscillate near (— vr, 0) and (0, vr). We, then, study the 
crossover between the two behaviors based on an analytical solution. 

One of the powerful methods analytically dealing with heteroclinic/homoclinic 
tangles of invariant manifolds is the method of asymptotics beyond all orders (ABAO 
method). The key idea is to employ the so-called inner equation and to investigate 
it with the aid of the resurgence theorjEl and the Borel resummation. The inner 
equation magnifies the behavior of the solution near its singularities and bridges the 
solutions in different sectors. 

The idea of the inner equation was first used by LazutkirP to derive the first 
crossing angle between the stable and unstable manifolds of Chirikov's standard map, 
and by Kruskal and Segur,'^ to study a singular perturbation problem of ordinary dif- 
ferential equations. Lazutkin and his 

coworkerPE»ee developed their method, and 
Lazutkin, Gelfreich and Svanidz^^ proposed a method of systematically improving 
the solution of the inner equation for the standard map. The first rigorous proof 
of the method was given by Gelfreich.EBCl He and his cohaboratorfPlC^'ESHIl 
have also studied the splitting of separatrices and related bifurcations for various 
systems. Meanwhile, Hakim and 

MalliclPZJ 

used the technique of the Borel transfor- 
mation within this context and obtained the first crossing angle previously obtained 
in Ref. EI). Tovbis, Tsuchiya and Jaff^^ improved their method and derived ana- 
lytical approximations of perturbed unstable/stable manifolds for the Henon map. 
Nakamura and HamadcP^'l^ studied the separatrix splitting for the cubic map and, 
as in the work of Ref. llUj) . Nakamura and KushibJ^ considered the systematic im- 
provement of the solution of the standard map. Also, ABAO method was applied to 
some higher dimensional systems by Hirata, Nozaki and Konishi.'^''^ 
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Roughly speaking, the procedure obtaining an approximate expression of the 
unstable manifold can be summarized as follows: 

(i) Find singularities in the complex time domain of the lowest order perturbative 
solution. 

(ii) Magnify the neighborhood of a singularity closest to the real axis and derive an 
asymptotic expansion of the lowest order solution valid in a certain sector. 

(iii) The Borel transformation and the resurgence theory are used to construct the 
asymptotic expansion which is valid in the other sector. Usually, there appear 
additional terms which are exponentially small for real time values. 

(iv) Going back to the original equation, derive equations for the exponentially 
small corrections corresponding to the terms found in the previous step. An 
appropriate solution is chosen by matching its asymptotic expansion with that 
obtained in the previous step. 

In this paper, with the aid of the ABAO method, the reconnection of the unstable 
manifold is studied for the Harper map. Our presentation mainly follows the work of 
Nakamura and Kushibe.l^ In the text, an approximate expression of the unstable 
manifold is constructed and only the result is given for the stable manifold. 

The rest of this paper is arranged as follows: In the next section, an outline of 
our result is presented without proofs and the details, which are given in Sees. 3-6. In 
Sec. 3, a perturbative solution is constructed by the Melnikov perturbation method 
and its asymptotic expansion is derived. In Sec. 4, the resurgence theory and the 
Borel resummation method are applied to solve the inner equation, and asymptotic 
expansions of the additional terms are derived. In Sec. 5, the exponentially small cor- 
rections are obtained by matching the solutions of the original and inner equations. 
As will be discussed later, the separatrix solution of ()l-2() has singularities along two 
lines both parallel to the imaginary axis and, among them two pairs on individual 
lines mainly contribute to the heteroclinic tangles of the manifolds. Sees. 4 and 5 
deal with one pair of singularities which have smaller real part and, in Sec. 6, the 
same analysis is repeated for the other pair of singularities. In Sec. 7, the derived so- 
lution is compared with numerical calculations and the reconnection of the unstable 
manifold is discussed. Sec. 8 is devoted to summary. 



In this section, we summarize the construction of the unstable manifolds of the 
Harper map with the aid of the method of asymptotics beyond all orders (ABAO 
method). The details of the analysis will be described in Sees. 13161 

2.1. Unstable manifold 

We are interested in the unstable manifold of the Harper map satisfying the 
boundary condition: 



When the parameter a is small enough, the solution of ()l-2() with the same boundary 
condition provides a good approximation for t <^ 0. However, the solution is not 
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valid for the whole real t since exponentially small terms should be added when t 
exceeds certain values t = Ti,T2 (Ti < T2). In fact, the unstable manifold is well- 
approximated by a double expansion with respect to a and e = a^e~'^ (c > 0,j: 
some constant): 

vit)\ _ fvo{t,a)\ ^.f. /Re(£;i(tr.«(t,a) 

where S{t) denotes the step function, Vn\t,a), Un\t,a), Vn\t,a,e) and Un\t,a,e) 
admit power series expansions with respect to a, and Ei{t) is an exponentially small 
but highly oscillatory function of t (more precisely, Ei{t) = 0(e) and is periodic with 
period a). 

The difference equations for Vn\un'^ are obtained by substituting (|2-2|) into (1.1) 
and equating the terms with the same powers in Ei{t). Particularly, vf\uf^ obey a 
linear equation and their values at t = are necessary for fixing them. Such bound- 
ary values are provided with the aid of the Borel summable asymptotic expansions 
valid near a singularities in the complex t-plane. To see the method more in detail, 
let us focus on the time domains — 00 < t < Ti and Ti < t < T2. 

2.2. Melnikov perturbation 

The solution fo(t),uo(t) describing the unstable manifold in the domain 
—00 < t < Ti is given by the conventional Melnikov perturbation: 

00 00 
VQ{t, cj) = ^ o"*i;oi(t), uo(t, 0-) = X] cF'uQi{t) 

i=Q i=0 

The lowest order solutions uoo(i), ■yoo(^) are nothing but the solutions of the differ- 
ential equation Hl-2() with the boundary condition H2-l() . After fixing the origin of t 
by vq{T) = 0((j^), the terms (uoo(*)) woo(O) and {vQi{t),uoi{t)) {i > 1) are found to 
have, respectively, logarithmic branch points and poles at 



(n = 0,±l,±2,---) (2-3) 



where T = In ]Jy^ ^ ^' other hand, by fixing the time origin as fo(0) = 

7r/2 + 0((T^), the perturbative solution of the unstable manifold for A; = 1 have 
singularities at 



r 

n -\ — \ iri 
2 
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As one can easily seen from the definition of T, T diverges as ^ 1 — 0. We choose 
the time origin so that the solution converges as ^ 1 — and the location of the 
singularities in the limit of A; ^ 1 — matches with that for the case of k = 1. We 
remark that the choice of the time origin, vo{T) = 0*\ gives the solution which does 
not converge as /c ^ 1 — (see appendix IE)) . 

In the domain Re t < of the complex time plane, (fo(t), uoit)) is a reasonable 
approximation of the unstable manifold since individual terms of their expansions 
have no singularities there. The expressions valid in the domain Re t > can be 
constructed by analytically continuing vq , uq along any curve which does not cross 
the lines t = zizis {s > -^^) emanating from the singularities of vq,uq. Then, it 
is convenient to study the analytic continuation along the curve passing near the 
singular points since the divergence of the perturbative solutions is expected to 
magnify the small corrections. 

2.3. Behavior near a singularity 

Now we focus on the behavior of the perturbative solution near a singular point 
t = ti = The corrections voi,uoi have the Laurent expansions: 



a 



^o^{t)) {t-t,Yf^\h?{t-try 



Hence, higher order corrections are more singular. In such a case, the behavior of 
the most divergent term is important. Because the most divergent contribution, e.g., 
to vq has the form:^^ (T*o[*V(* ~ it is convenient to introduce a new variable 
z = and to rearrange the expansions as double expansions with respect to a 
and z. Since one has 

(T^voi{t,a) = ^ +CT^ + --- 
z z" 

z^ z^ 
a^vo^{t,a) = ^ + a^ + ■■■ 
and the similar expression for uqi, the perturbative solution can be reexpressed as 
Voit,cr)\ i f V0iit)\ _ f Vooit)\ \^ I fVoiiz)\ 



where Vqi , Uqi are given by 

c 

^oK^)=EiT' Uoiiz)=^^ (2-6) 



. , Z -. z 

l=\ 1=1 



*■* This choice gives the symmetric unstable manifold with respect to u axis. But it makes Stokes 
multiplier divergent and it indicates that the unstable manifold is not symmetric with respect to it 
axis even within the perturabation approach. 
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We remark that these series are divergent and are meaningful only as asymptotic 
expansions. 

In short, the perturbative solution valid in the domain — oo < Re t < turns out 
to be an asymptotic expansion near its singularity. Therefore, a careful treatment is 
necessary to analytically continue the solution across the line Re t = 0. 

2.4. Borel resummation and analytic continuation 

Fortunately, the series (|2-6|) are Borel summable. This can be shown as follows: 

(i) The equations for the Borel transforms Vqi (p) , Uqi (p) oiVoi{z), Uqi (z) are derived 
from the original equation and the expansion H2-5|) . 

(ii) The singularities of the Borel transforms are found with the aid of the linearized 
inner equation. Indeed, according to the resurgence theory,'^'^ the functional 
forms of the Borel transforms Vqi (p) , Uqi (p) are obtained from the linearized 
equation except some constants. 

(iii) The Borel transformed equation is solved. This is carried out as follows: The 
previous step determines the functional forms of Vbz(p), Uqi{jp) up to a few ar- 
bitrary constants. These constants are settled by comparing the power series 
expansions of the guessed solutions with the power series solutions of the Borel 
transformed equation, which are numerically derived. 

Then, the Borel transforms Vqi{p)^Uqi{p) are found to have singularities along the 
lines p = is, {s < — 27r or s > 2tt) and are analytic near the origin. The analyticity 
implies the existence of the asymptotic expansions of Voi{z), Uqi{z) at |z| ^ oo in a 
certain sector. Hence, the divergent series H2-6() are Borel summable. 

Now we consider the analytic continuations of Vqi{z),Uqi{z). Since negative 
powers of z in Vqi, UqI are expressed as the inverse Borel transformation (namely 
the Laplace transformation) of functions Vqi (p) , Uqi {p) which are analytic except on 
the lines p = is, (s < —2tt or s > 2-k), their analytic continuations from the sectors 
Re z < 0, Im z < to Rez > 0, Imz < 0, are obtained by appropriately changing the 
integration path of the Borel transformation. Then, in the sector < Re z, Im z < 0, 
the solutions are expressed as 




where Vii{z),Uii{z) are additional contributions and only the terms dominant for 
real t are retained. More concretely, up to and one has 




where 

Aa = A^ + aA'^ + --- {a = A, B) 
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and the coefficients Af are numerically determined. Note that the analytic contin- 
uation 1)2 -Tf) simply provides the asymptotic expansion near t = ti and that further 
investigation is necessary to construct the analytically continued solution. 

2.5. Solutions in different domains 

Now we consider the solution for Re t > 0. By comparing (|2-2j) with (|2-8|) and 
reminding az = t — ti, one has 

Ei{t) = e-— , (2-9) 

a 

and, thus, e = e /a. Because Ei{t) is exponentially small with respect to a and 
periodic in t with period cr, equations for {vn\un'^) are obtained by substituting 
(|2-2j) to and equating the terms with the same powers in Ei{t). Substitut- 

ing Vn\t,cF) = J2i^o ^""^ni i't) ^^^^ Un\t,a) = X]^o (*) ^lie so-obtained 
difference equations, one has a system of differential equations (n, i) ^ (0,0): 

-i'^'W +E (I)' %^ = --iVW-snoo(t) + 

^7^ • (1) (2-10) 

where and are functions of Vj\} {j < n, k < i), v^^j! {k < i), vok {k < i) 

and uj'^ {j < n, k < i), u~^l {k < i), uok {k < i), respectively. The general solution 
of the system involves several constants of integration, which should be chosen so 
that the matching condition: 




(2-11) 



is satisfied for small a and large z. Another matching condition at t = t\ similar to 
(|2-lip is obtained by repeating the procedures explained in the previous subsections. 
One can prove, from real analyticity of (woj, UQj), that the additional term obtained 
from the anlysis of t = is a conjugate to the term obtained from the analysis 
of i = ti (see appendix . These two matching conditions fix the constants of 
integration and, thus, the solutions Vn\t,a),Un\t,(j). As easily seen, the solutions 
have no singularity in < Re t < 2T. 

This completes the derivation of the solution valid in the domain < Re t < 2T. 
When Re t exceeds 2T, new terms will appear and the solutions valid in the 
domain Re t > 2T can be obtained as follows: 
(i) The asymptotic expansions valid near t = t2 = 2T + derived as a 

function of z' = {t — t2)/(T and are analytically continued from the domain 
Rez' < 0, Imz' < to Rez' > 0, lm.z' < with the aid of the Borel resummation 
and the resurgence theory. 
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(ii) The system of differential equations for the additional terms in Re t > 2T is 
derived from (|l-ip and ()2-2() . And its general solutions are calculated. 

(iii) From the matching conditions such as H2-ll() near the singularities t = ^2)^2' 
the constants of integration are fixed and one gets the desired solutions. 

The final solution is shown to be analytic in the domain Re t > 2T and, thus, the 
description of the unstable manifold on the whole time domain is completed. Note 
that, when t passes the line Re t = 2T, the terms added in the domain < Re t < 2T 
generate new terms, which are of order of and negligible for real t. It will be proved 
in Sec. 6. 

2.6. Approximate unstable manifold 

Up to the order of and cr^e, the unstable manifold is expressed as 



Vu{t)= Vooit) + cr^U02(t) 



S{t)Re 



4yl(l) 



-e " 



~X2{t)e~ 



+ S{t - 2T)Re 



e— { f ^ TT^, 1 xi{t) + X2{t) \ e 



Ak{k - 1 

Uu{t)= uooit) + + a\o2{t) + Qn'o2(t) - ^l/i 

+ S{t)Re 



— — e— y2{t) + ) e"" 



+ S{t - 2T)Re 



la V 2 

4yl(2) 2.u^ f /T(fc- 1)2 - (1 + fc) 



Ak{k - 1) 



yiit) + a 



y'lit) 



+ {y2it) + a 



y'2it) 



_ 2-Kit 

e " 



where S{t) stands for the step function and the other functions are defined as follows: 



'^'oo(^) = ~2tan ^ sinh — VA; cosh 

Vk 



itoo(i) = 2 tan 



-1 



vm{t) = -- 



sinh ^/kt + cosh y/kt 

( sinh y/kt — \fk cosh \fkt ) ( cosh \fkt — y/k sinh \/l:t 
c[{t)+xi{t) { kt-2Vk^ 



(2-12) 



U02{t) 



24 



2y[{t) - yi{t) { kt - 2Vk- 



(1 + k) cosh^ ^/kt — 2\fk sinh \fkt cosh \fkt 
\ sinh \fkt — \fk cosh \fkt ) ( cosh \fkt — \fk sinh \fkt 



(1 + /c) cosh^ \fkt — 2y/k sinh y/kt cosh 



(2-13) 



xi{t) = -2Vk 
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1 + ^sinh y/kt — y/k cosh 



, , sinh \/kt — \/k cosh vfct 
yi{t) = -2fc z -2 (2-14) 



X2it) 



1 + ^sinh — cosh j 

- — ('sinh Vkt — \/k cosh ^fkt^ H — 2 

^ ^ 4 V/c 1 ^ (^gij^j^ Vfct - Vfc cosh V^i) 

1 — k f iir \ , , 
2/2 (i) = ^ ( cosh Vkt — sinh H — 7-7^ ; — 2 



4^ ^ ^ 4^/fc 1 ^ J^g-^j^ ^t-^/k cosh V^t) ' 



1 — /c / iiT \ , , 



The coefficients A^'^\ are numerically determined as 

= i47r3^i + avr^ I - (fc - 1)^2 - ^^^^4 
yl(2) = _i47r3^i + CTvr^ |-(A: - 1)^2 + ^^^^2^^^ ''^ '^^^ 

where Ai = 0.27893, B2 = 0.14, = 3.503. 

§3. Melnikov Perturbation 

In this section, we start with the investigation of the perturbative solution, 
which provides a good approximation of the unstable manifold for t — > — 00. The 
lowest order terms in the expansion VQ{t^ a) = X^^q '7*^0i(i)) ^^0(^5 o") = Yl'^o (^^uoi{t) 
satisfy Hl-2p and the first and second order terms obey 

Uoi(t) + uoi{t) cosuoo(t) = - ^VQQ{t) 

u'oiit) - kvoi{t) cos voo{t) = ^noo(t) (3-1) 

VQ2(t) + uo2(t)cosuoo(t) = -^VQi{t) - -VQQ{t) + — sinuoo(t) 

1 , 1 (t) 

«02(*) - f^V02{t) COS voo{t) = -Moi(t) - -noo(t) - k sini)oo(t) (3-2) 

The lowest order solution of the unstable manifold with ^00(7") = is given by 
H2-12p . Because of the boundary conditions fon(i) — > 0, UQn{t) for i — > — 00 and 



(2-15) 



10 



S. Ajisaka and S. Tasaki 



t = —oo 




Ret 

~t = ao 



X : pole 
— :branchcut 

Fig. 2. Singular points (crosses) of the lowest order solution vq, uq in the complex time domain. 



voi{T) = vo3{T) = 0, ^02(7") = ^Yg-Vl — k, the first and the third order solutions 
are found to be 



voi{t) = 0, uoi{t) = -yi{t) 



vo3it) = , uo3{t) = ^n'o2(t) - j^y'lit) 

and the second order solution is given by (|2-13p . It will be discussed in appendix lEl 
we fix a time origin so that the final solution has no divergent term and the choice 
forbids a symmetric perturbative solution with respect to u axis. We observe that 
each order solution has two sequences of singular points in the complex t-plane (see 
Fig.EJ. 



^ f n + 



0,±1,±2, 



(3-3) 



where T = In > 0. This implies that the perturbative solution approximates 
the unstable manifold well in the domain: Re t < 0. 

In the rest of this section, we examine the behavior of the perturbative solution 
near the singular point in the upper half plane: ti = closest to the real axis. By 
substituting the Laurent expansion of each perturbative term to the set of equations 
(|2-10j) . we inductively get 



a a, 



(i) 



+ 



a a 



i-l 



+ 



(i>l) 



and, thus. 



vo{t, a) - voo{t) = ^ a'voiit) = X] X] 



i=l 



i=l 1=0 



a 



(3-4) 
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This indicates that higher order terms with respect to a have higher order poles at 
t = ti and that a careful investigation is necessary. For this purpose, it is convenient 
to introduce a scaled variable z = {t - ti)/a^^^^^ Then, one finds that 
Vo{z, a) = vo{ti + az, a) - voo{ti + ctz) and Uo{z, a) = uo(ti + az, a) - noo(ii + (tz) 
can be expanded into power series with respect to a: 



Vo{z,a) 
Uoiz,a) 



Sri fVoliz) 
^ ' Uoiiz) 



1=0 



with 



Voi{z) 
Uoiiz) 



oo ^ 



(3-5) 



i=l 



and that (Voo,C^oo)) (^oi)t^oi)i ••• correspond, respectively, to the most divergent 
term, the next divergent term ... at t = ti. It turns out that the expansions H3-5|) 
of Vol , Uoi are Borel summable asymptotic expansions which are valid in the sector 
Re z < 0. Hence, when one crosses the line Re z = 0, new terms would be added. 
These aspects will be discussed in the next section. 

§4. Inner Equation and Analytic Continuation of Solutions 



4.1. Inner equation 

In terms of the variables V{z,a) = v{ti + az,a) — foo(*i + i^z) and 
U {z, a) = u(ti + az, a) — uoo{ti + az), the Harper map (|1-1(1 reads as 

AV{z, a) = -a sin {U {z, a) + uoo(*i + (^z)) - Avooih + crz) 
AU{z -l,a) = kasm{V{z,a) + voo{ti + crz)) - Auoo{ti + a{z - 1)) (4-1) 

where A stands for the difference operator: Af{z) = f{z + 1) — f{z). This equation 
is used to study the asymptotic property of the solutions for \z\ — > oo *^ and will be 
referred to as an inner equation. On the other hand, the map l|l-l|) in the original 
time variable is referred to as an outer equation. 

Let Vo{z, a), Uo{z, a) be the solution of ()4-l() corresponding to the perturbative 
solution and let 

'Vo{z,a)^ ^^^ifVoiiz)' 
Uo{z,a) 



1=0 



\Uoiiz) 



(4-2) 



be their expansions with respect to a, then the equations for Vqi, Uqi are obtained 
by substituting H4-2|) into 1)4- 1() . The first two sets are 



AVooiz) 
AUooiz) 
AVoiiz) 
AUoiiz) 



-iUooiz) 

z 



Z+l 

■ k-l , Uoiiz) 



iln[ 1 + - 

z 



+ iln[ 1 + - 

z 



z 



^Ol(^ + l) 
Z+l 



1-k 



,iyoo(.+i) + ^(l_^) 



(4-3) 



(4-4) 



*■* The asymptotic behavior for \z\ — > oo corresponds to that for t ^ fi in the original varibale. 
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From these equations, one can directly obtain the expansion (|3-5|) discussed in the 
previous section. Indeed, by posing the matching condition at ^ = 0: 



Res 



\Vok{z) 



= Res 

z=0 



Vok+lit) 
Uok+lit) 



(4-5) 



t=ti 



one has 



^°°^")= 12;^ -4320i^ + °b 

i i i 191i / 1 

^oo(^) +^ +^ + O 

... X i A:ti + 1 I n 1 



- 1) , i + 1) i k{ti + l) I 
^''^'^ = +24^ 48^^ +° i ^ 



In the following subsections, we shall show that these expansions are Borel summable 
and that their analytic continuations from the sector Re z < 0, Im z < to 
Re z > 0, Im 2; < give necessary information for the construction of the ap- 
proximate unstable manifold for < Re t < 2T. Note that , one can prove from the 
above asymptotic expansions and (|4-3|) that all the coefficient of 1/z exponents of 
Voo(-z), Uqq{z) are purely imaginary. 

4.2. Borel transformation and analytic continuation 
We define their Borel transforms Vo(j'), Uq{p) by 



Foo(^) = L[Vo{p)\{z) ^ \ dp e-P' Vo{p) 

Jo 

Uoo{z) = L[Uo{p)]{z) = / dp e-P' U^{p) (4-6) 



The integration path is so chosen that the left hand sides are regular for 
Re z — > —00 since we are investigating the unstable manifold. As easily seen from 
631), Vq{p), Uo{p) satisfy 



i{e-P-l)Vo{p) = l+ dp\Y, 



{-^Tui''''\p)\ 1 



1^1 ™' J P 



i(l - e^)mp) = 1+ Tdplf^ , - ^ (4-7) 



nl J p 

where V^*^^ denotes the nth convolution recursively defined by 

Vi*''\p)= rdxVo{p-x)vi*^''^''>\x) (4-8) 



Reconnection of Unstable/Stable Manifolds of the Harper Map 



13 



cut 




x:pole 



cut 



Fig. 3. Integration path in p-plane. 

The factors (e~^ — 1) and (1 — e^) of ()4-7|) indicate that the Borel transforms 
Vq{p), Uq{p) may have singularities at p = ib27rnz (n = 1,2, • • • ) and branch cuts 
starting from them (cf. Fig. 

Then, the analytic continuations of Voo{z) and Uooiz) from Re z < 0, Im z < 
to Re z > 0, Im z < are obtained by rotating the integration path in the p-plane 
from C to C as shown in Fig. |31 Note that Voo{z) and Uoo{z) should be analytically 
continued along a curve passing below the singularity z = where they have no 
cut (cf. Fig. 121). This corresponds to counterclockwise rotation of z, e.g., along 
a semicircle below z = and, thus, one should rotate the p-integration path in a 
clockwise way (cf. Fig. In short, the analytic continuation is given by 



^00 (^) 



dpe-P'Vo{p) 



dp e-P' Voip) - / dp e-P' Voip) (4-9) 



where 6 = tt/2 + e' with e' being a small positive number. The first term has the 
same asymptotic expansion with respect to z as in the domain Re z < 0, Im z < 
and the second term is the desired additional term. 

4.3. Linearized equation 

If one can explicitly write down the Borel-transformed solution of the inner 
equation, one has only to change the integration path. But it is difficult to find the 
location and order of singularities in a straightforward way. Instead, we investigate 
them with the aid of the linearized inner equation. 

As seen from the Borel transformed inner equation H4-7() . singularities of the 
Borel transforms of Voo{z) and Uo(){z)neavest to the real axis are located at p = ib27rz. 
According to the resurgence theory ,1^1 'El the order of these singularities is given by 
the solution of the linearized equation of ()4-3() : 



A^{z) 



-iUooiz) 



A^{z) 



-^{z + l) . 



(4-10) 



z z+1 
More precisely, let {Va{z) , Ua{z)) and {Vb{z),Ub{z)) be two linearly independent 
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solutions of (|4-1U() and let {Va{p), {7a (p)) and {Vb(j>), Ub{p)) be the functions satis- 
fying 



7 

A „^ J aB 



then, with appropriately chosen constants Aq and Aq , the Borel transforms of Vqo{z) 
and Uoo{z) are given by 



(Voip)\ 



^A (Va{p)\ , ^B (Vb{p)\ , (Rvip) ) .4.12 

^° [ua{p))^ ' + ^ 



in a neighborhood of a singular point p = liri, where the functions Rv{p) , Ruip) are 
regular near p = 2Tri. 

With the aid of this observation, the Borel transforms Vq{p) and Uq{p) are 
determined as follows: At first, with the aid of the formulas*^ 



dp e-P^ 
dp e-P' 



U - 1)! 



7 



{p — 2TTni)^ 

(p — 27riny \n(p — 27rm) 



27ri(-z)J'-ie-2™^ 
-27rz 



and the relation (|4-11|) . one derives approximate functional forms of Va{p), Ua{p) 
(a = A, B) near p = 27ri from the asymptotic 1/z-expansion of the solution of ()4-10() . 

As pointed out at the end of Sec. 14.11 all the coefficients of 1/z-expansions of 
Vqq{z) and Uqq{z) are purely imaginary and this implies that Vb(p), Uq{p) are purely 
imaginary for real p. By taking this property into account, the functional forms of 
^o(p)) Uq{p) are guessed from H4-12() . Next, the power series solution of the Borel 
transformed inner equation 1)4 -7^ is obtained numerically. Finally, the coefficients 
Aq and A^ are determined by comparing the power series solution of Vq{p), Uq{p) 
and the power series expansion of their guessed functional forms. This calculation 
will be carried out in the next subsection. 

In the rest of this subsection, we give a plausible argument why the Borel trans- 
formed inner solution near the singularity can be determined from the linearized 
equation. We restrict ourselves to the case where the Borel transformed solution has 
the form**): 

= Y^ ( 9niP - 27rni) + hl{p - 27rni) ln{p - 27rm)\ 
^^ygilip — 2'Kni) + h'^{p — 27rni)ln{p — 2TTni) J 

where g^ip) {a = u,v) are meromorphic with a pole only at the origin and h'^{p) 
{a = u,v) are analytic and e~^^/i"(p)(Im z < 0) are rapidly decreasing for p — > ioo. 

*■* We have chosen the branch cut of the logarithm along the positive imaginary axis as shown 
in Fig. 121 Note that this choice fixes the Stokes line on the negative imaginary axis in the z-plane, 
which corresponds to the line in the original t-plane joining ti and its complex conjugate. 

**' For the Henon map, Gelfreich and SauziiP have shown that the solution of the Borel trans- 
formed inner equation has such a form in a neighborhood of the singularity. 
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Then, their contour integrals along the path 7 are given by 



dpe 



-pz 



(Vo{p)\ 



n=l 



-2ninz 



where 



(Vno{z)\ ^ 27r 



e^^^^^Res 



-pz 



+ 



(Vnoiz) 
\Un0{z) 



dqe 



(4-15) 



p=0 



-iqz 



(.4-16) 



The first terms of the right hand side of (|4-15() are polynomials of z and the second 
terms can be expanded into a power series with respect to 1/z. On the other hand, 
since 



Voo{z)- f^dpe-P^VoipY 
Uoo{z)-f^dpe-P^Uo(j))^ 



Vooiz) 
Uooiz) 



n=l 



-2ninz 



Vnoiz) 
Uno{z) 



and (Voo(-2), Uqq^z)) are solutions of ()4-3() . by substituting the above expansion into 
(|4-3() and comparing the terms with the same powers in e~^'^*^, one finds that Vio{z), 
Ui(){z) satisfy the linearized equation (|4-lUp and, thus, are linear combinations of 
{Va{z),Ua{z)) and {Vb{z),Ub{z)). As a consequence, there exist constants Aq and 



Aq such that 



J 7 V' 



glip - 27rz) + h\{p- liri) \n{p - 2m) 
gf{p - 2TTi) +h{{p- 2m) ln(p - 27ri) 



(Vio{z) 
\Uidz) 



-2niz 



^^0 



VAiz) 
Ua{z) 



B 



Vb{z) 
Ub{z) 



(4-17) 



This relation corresponds to 1)4 -1 2(1 and implies that the functional form of the Borel 
transformed solution Vq{p), Uq{p) can be determined by the linearized equation. 

4.4. Determination of Vo{p), Uo{p) 

Now we calculate Vo(p), Uo(j>) following the procedure described in the previ- 
ous subsection. Firstly, the linearized equation (|4-11() has two linearly independent 
solutions {Va{z), Ua{z)) and {Vb{z), Ub{z)), which can be expanded for large z as 



VAiz) = - 
Vb{z) 



1 1 

z bz"^ 

1 

= z + 

540z3 



^-(^) = ^-2;^ + T2;^ + -- 

Ub(z) = z-\ \ K 

^ ' 2 \2z 24^2 



+ 



11 



2160z^ 



+ 



(4-18) 



Then, (|¥T^ leads to 



VAiv) 
Vb{p) 



2TTi 

1 



ln{p — 2'ni) + ■ 



1 



27ri {p — 27ri)" 



+ • 



(4-19) 
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where the omitted terms correspond to z ™ (m > 2). Hence, one has 

where the rest part is either regular or less singular at p = 27ri. As mentioned at 
the end of Sec. 14.11 all the coefficients of the l/z-expansion of Vbo(z) are purely 
imaginary and, thus, Vq{p) is purely imaginary for real p. This property determines 
a part of the regular components of Vo{p) and, as a result, one can guess 



1 [ -^o" , < 

Im I (p — 27ri)2 (p + 27ri 

Re< /(^)(p) + Im< - Re.l^/l?(p) - W^/i\)(p) + • • • (4-21) 

where the following auxiliary functions are introduced 



fj""^ (P) - ( _ 2^i)i+i + + 27ri)J+i ) - 0' 1' 2' • • • ) 
f('\p) = 2L( I I ^ (. = 1 2 •••) 

= ^ [ln(p - 27Ti) + ln(p + 27ri)] 

/i'l (P) = -7^ - 27ri) - ln(p + 27ri)] 
zvr 

We note that the above auxiliary functions have power series of p with purely imag- 
inary coefficients. One can easily show that the neglected terms of (|4-21|) are either 
regular or less singular at p = zizlni. As easily seen, (|4-21[) admits the power series 
expansion near p = 0: 

n=l ^ 

where M = — ^^™^o _j_ ^^^Jp — log(47r^). In a similar way, one can guess 

UM = M' + V I (2„ + 2) - '-(f + ^°7'^) + (-1).. ( A)^-' 

uv/^^ ^1 47r3 ^ ' 7r(2n+l) iir^ ( ^ ' \2ttJ 



n=0 

oo 



where M' = M + (^^;^+<A^) _ 
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Now we numerically solve the Borel transformed inner equation (|4-7|) . By sub- 
stituting the power series expansions: 



oo 



^o(p) = XI ""^P"^' Uq{p) = X 6np"; ao = 0, 60 = 

n=0 n=0 

into ()4-7|) and comparing term by term, the coefficients 0^, 6„ are determined recur- 
sively and one has 

za2„+i(-l)"(2vr)2"+i ^ -(2n + 2)^i + ^ 00) 

2n -|- 1 

ia2n = 0, (Vn) (4-24) 

z62n+i(-l)"(27r)2"+i ^ -(2n + 2)^i + -if-r ^ 00) 

2n -|- 1 

i62n(-ir(27r)2" ^ ^3 (n^oo) (4-25) 

where Ai = 0.27893, A2 = 0.417, ^3 = 0.87628. 

By comparing with and with one finds that Aq and 

yl^ are purely imaginary and 

< = 47r3iAi , + ^ = -^^ = vriAa , < = 47r2z^3 . (4-26) 

Because of the first and third equalities of ()4-26p , the ratio ^3 /Ai has to be t^Aq / Aq = 
IT and the evaluated values give an excellent agreement: A^/Ai = 3.1416. Also, the 
first and second equalities require {QA2)/{tt'^Ai) to be unity and the evaluated values 
are consistent: (6A2)/(7r2^i) = 0.909*). 

Here we remark about the contributions from a singularity at p = ipc {Pc > 27r). 
If the singularity is of the same type as that at p = 27ri, the contributions to the 
coefficient of are exponentially smaller by a factor of (^|^^ than those from 
p = 2TTi. Thus, such contributions do not affect the present analysis and, at the 
same time, are difficult to be evaluated. 

In short, following the procedure discussed in the previous subsection, we have 
obtained the Borel transforms Vo{p), Uo{p): 

Mp) = -^f^Il{p) + ^f['\p) 

Uoip) = -^f^Jhp) + ^ [fi'\p) - '^fPiP)) (4-27) 



*-* The numerical evaluations of Ai, A3 are quite robust. On the other hand, A2 is sensitive 
to an error of Ai since the coefficients of l/(2n + 1) are fitted after subtracting the leading terms 
of order n. However, we observe that the coefficients of l/(2n + 1) in a2„+i( — l)"(27r)^"+^ and 
62ti+i( — l)"(27r)^"^^ are always the same, or the second equation of (14-2611 holds. And the value of 
A2 are fitted so that the first and second equations are well satisfied. 
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4.5. Borel transforms o/Voi(z) and Uqi{z) 

Let V^{z)^ Uqi{z) be the sum of negative powers of z in the asymptotic expan- 
sion Vqi{z), Uqi{z), and let Vi{p)^Ui{p) be the Borel transforms of Vq^{z), Uqi{z) 
in the sector including z = — oo, then they provide the following additional terms in 
the other sector including z = +00 

dp e-P^ ( '-'^f'l] = e"2-- f;^!!^^!^ + • • • (4-28) 




Uii{z) 

where the neglected terms are proportional to a factor of e~^^''^ {pc > 2tt). 

Similarly to the construction of Vo(p), Uo{p), one can get Vi{p), Ui{p) from 
Vii{z), Uii{z), which obey: 

z z 
AUuiz - 1) = Mf)±^lfMf),^v^oo(.) + /io(^)Joi(^)^.Voo(.) (4.23) 

where Vio, Uiq are introduced in (|4-15p . 

Because of the linearity of (|4-29jl . Vn and Un are given by 

Vuiz)\_ A (V,iiz)\ B (Vf,{z)\ ^ fVA{z)\ B (Vb{z)\ 

where the constants Aq (a = A, B) are introduced in the previous subsection, yl" 
[a = A, B) are new constants, and V-yi{z), Uyy{z) (a = A, B) are solutions of 

z z 
AUt,{z - 1) = ]^M±^!f]:i^e^v.oo(.) + /»(^)Joi(^),.Voo(.) (4.31) 

such that the asymptotic expansions of ^^1(2) (a = ^, -B) have no terms proportional 
to z nor 1/z. Then, asymptotic expansions of V\\(z\ \J\\(£) are obtained by exactly 
the same procedure as in the previous subsection and, for large \z\^ we have 

where the terms of order z~™ (m > 0) are omitted since the numerical estimation 
of their contributions to Vi(p)) Ui{p) is very difficult*^. 

Now we go back to the equation H4-4() of Vbi(z), Uoi{z). Because it is linear, one 
gets (See appendix O foi^ more details.) 

V,-,iz)\ _ ^ / kt, + l (Va{z)\ ^ k-l / Vb{z) \\ (4.33) 



UqM) I 24 \Ua{z)) 2 \Ub{z)-{z-\)^ 

*■* Such an estimation requires a separation of the terms of order 1/n from the coefficients a„, 
fen of the power series expansion of Vi,Ui. But, since their leading order terms are proportional to 
, it is very difficult. 
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As one can see from the definition of | J^l^) I, the Borel transform of it has real 

\Ua{z)J 

values for real p. On the other hand, H4-3'2[) shows that Vi{p), Ui{p) have the first 
and second order poles at p = 27ri. These observations imply 

(IS) ^ - '^'A-^) - i^^'"' ' " - '"'^^) (!) 

-nHk-l)B,fl'\p)(J'j (4-34) 



j;;;g).4.',*-.)B,-(_\ 

The constants Bi, B2, -B3 and -B4 can be numerically evaluated by considering 
the Borel transformations of Va, Ua, Vb — z and Ub — z — 1/2 (for more detail, see 
Appendix^. Then, we have 

Bi = 0.01480, B2 = 0.14, B3 = 0.186, B4 = 3.503 (4-36) 

By substituting ()4-35|) into H4-32() . one gets 

= 2ATT^iBi = Gn^iBs 

Af = -(k - l)^^B, - ^^7r=^54. (4-37) 

Eq. H4-37() and Aq = i47r^^i imply two relations among Ai, B2 and B4, which are 
satisfied by the present numerical estimations rather well: 

Ai 2Ai 

0.99987 ~ 1 — = 0.99974 ~ 1 . 



6irBi 37ri33 
Hence, we obtain 

m=^^^^^fi'\p)+Aff[^\p) 

UAp) = {f^'\p) - f['\p)) + Aff[^\p) (4.38) 

4.6. Asymptotic expansions of the additional terms 

In summary, we have shown, up to the order of a^, that V{z, a), U{z, a) acquire 
new terms in the sector including z = +00 

V{z,a)\ ^ ( Voo{z) + aVoi{z)\ 2.^z ( Vio{z) + aVn(z 



Uiz,a)J - \Uoo{z)+aUoi{z)J^^ ^'yUMz) + aUn{z))^'" ^^'^^^ 
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where the first term corresponds to the perturbative solution studied in the previous 
section and the second term is the additional term: 



Vio{z) + aVu{z) ^ ( + A^z\ +a^Afz + Ao — 



Ulo{z) + aUn{z) 



z +"^0 r+2 + 12z 



+ (7 



A?z-A^- ^'-'^ 



(4-40) 

{z' + z) 

(4-41) 



Although the functional forms of Flo, f^iOi ^ii) ^ii can be derived from their linear 
equations, the resurgence theory and Borel resummation play an essential role for 
the determination of the coefficients A^^A^^Af. 

In the next section, only the terms of order of z and are used to derive the 
solutions of the outer equation in the domain Re t > since cr/z-terms are not 
determined because of the nemerical difficulty. However, it should be noted that the 
existence of 1/z terms does show the fact that the point p = 27ri is a branch point of 
the Borel transforms of V{z, cr), U{z, a). Before closing this section, we remark that 
(|4-40|) can be rewritten as 



Vio{z) + aVii{z) 



-^ + {A^+aAf) 



z + a- 



k - 1 



§5. Matching between Inner and Outer Solutions 



5.1. Matching at the singular point t = ti 

In this section, solutions of the outer equations are constructed and they are 
matched with the inner solutions. We first consider the contribution from the sin- 
gularity t = ti. Corresponding to the expansion of the analytically continued inner 
solutions: V = Vq + e~'^^^^Vi + ■ ■ ■ , U = Uq + e~^''*^C/i + • • • , the original solutions 
V and u acquire new terms in a sector Re t > Re ti = of the ti-neighborhood: 

v{t) = vo{t, a) + vi{t, cr)e"^* H 

u{t) = no(t,cj) + Mi(t,(j)e"^* + --- (5-1) 

2 

_9 • 2TTi , 

Because of e = e e <^ , vi and ui are exponentially small with respect 
to a. By substituting (|5-H) into and comparing term by term, we obtain the 

equations for vi and ui: 

vi{t + a) — vi{t) = —aui(t) cos uo{t) 

ui{t + a) — ui{t) = kaviit + a) cos VQ{t + a) (5'2) 
Its solution is uniquely determined by the matching condition: 



_2Tni f vi(t, a)\ 
^ " \ui{t,a)J 
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or equivalently 



vi{ti + az,a) 
ui{ti + az,a) 



HiiHi f Vio{z) + aVu{z) + • 
\Uio{z) + aUniz) + 



This condition suggests the fohowing expansions: 

oo 

vi{t, a) = a^e^~ } vin{t)a'^ , ui{t, a) 



n=0 



a-'e <y 



oo 
n=0 



(5-4) 



ni„(t)cT" (5-5) 



where the exponent j will be determined later. Then, equations for v\q^u\q^v\\^u\\ 
are given by 



^io(0 = --"10(0 cosiioo(t) 
tiio(i) = A;fio(t)cost'oo(i) 



-uii (t) cos lioo (t) + liio (t)'uoi (i) sin uoo (t) 



- o^^ioW = ^ (^ii(*) cos^;oo(i) - t'io(i)^^oi(i) sinuooW) 



First, we solve H5-6|) . Since it is linear, one has 



fio(i) 
nio(t) 



+ 6 



X2(t) 



(5-6) 
(5-7) 
(5-8) 



where X2(t), ?/2(t) are defined by (PTl)) and lfTT3)) . 

On the other hand, up to the order of fj*^, the matching condition leads to 



vvd{az + ti) 
uio(cTz + tx) 



a- 



x\(gz + t\ 
yi{cFz + ti 

(Vio{z)+aVii{z) + 
\Uio{z) + cjUiiiz) + 



+ h 



X2{(7Z + ti) 

y2{az + ti) 



A, 



B Z 



dom. 



(5-9) 



where the subscript indicates to take the terms of order and 'dom.' stands for 
the largest part for z ^ oo. Because the left hand side of (|5-9|) starts from z, one 
should have j = —1. And the coefficients a and h are determined by the requirement 
that (|5-9|) admits well defined o" ^ limit: 







Therefore, up to the 0th order in a, we have 



Vlit) 

ui{t) 



2A^ 
( 

ia 



2Ai 



'(t-ti) / X2it) 



Similarly, up to the ffi'st order in a, one obtains 



ui{t) 



e <^ 

ia 



X2{t) 

y2{t)+ay'^{t)/2 



(5-10) 



(5-11) 



where A^^^ = Aq + aAf . Matching of higher order terms with respect to a is 
discussed in Appendix IfI 
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5.2. Contribution from the singular point at t = t\ 

So far, the behavior of the solution near t = ti has been considered. In order to 
derive the solution valid in the domain Re t > 0, the contribution from the complex 
conjugate singularity t = t\ should be taken into account. Because vo{t),UQ{t) is 
real analytic, it is simply the complex conjugate of the contribution from ti (See 
appendix El- Hence, the unstable manifold for t < 2T is described by 



Vuit) 
Uu{t) 



4yl(l) _2ni{t^ ( X2{t) 

y2{t)+ay'^{t)/2 



+ 5(t)Re 



where S{t) denotes a step function. 



§6. Contributions from Other Singularities 



(5-12) 



So far, contributions from ti and t\ have been considered. Here, we study 
contributions from other singular points and higher order terms. As easily seen, in 
the domain Re t > 0, the expression (|5-12j) of the unstable manifold has a singularity 

at 

Therefore, the solution would change its form when t exceeds 2T. As in the previous 
analysis, the additional terms can be found from the asymptotic behavior of the 
solution near t2 and its complex conjugate t^. 

The solution (|5-12|) valid vat < 2T is a sum of the perturbative solution and the 
terms proportional to and each term produces a new additional term in the 

sector Re t > 2T. We analyze these two terms separately. 
(1) The additional terms arising from the perturbative solution. 

The additional terms arising from the perturbative solution can be obtained 
exactly in the same way as those from t = t\ and we have 



f^o"i(^) 



kt2 



1 
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Va{z) 
-Ua{z] 



+ 



k-l 



-{Vb{z)- 
Ub{z) - {z 



(6-1) 



,(2) 



it) 



Re 



e <^ 



r(fc- 1)2 - (l + fc) 
Ak{k - 1) 

X2it) 

y2{t) + ay'2{t)/2 



Xl{t) 

yi{t) + ay[{t)/2 



(6-2) 



where yl(2) is given by yl*^2) = —iA-ji'^Ai + cttt'^ | — 



1)^2 + 



12 



(2) The additional terms arising from the term proportional to S{t). 



Let us prove that the terms added in the domain < Ret < 2T generate 
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new terms which are of order of and, thus, are neghgible. With the same 
analysis as for the most divergent terms from the perturbative solution, their 
most divergent terms near t = t2 are found to be 



ir 



a 



m 



VaIz') 
-Ua{z') 



iFa 2-Kt ( Vb(z' 

cos — ' 

2 a 



Ub{z')) ' 



(6-3) 



where i{t)^r are defined by 



t2{l - kf - 2(1 + k) 2-Kt ti(k-l) 2'Kt 

i{t) = — -4i \ ^ cos + ^ , ^ sin 

' 8k{k-l) a 8ik a 



r 



Vka 



(6-4) 



This can be represented with the aid of the Borel transformations: 

ir 



a 



-m 



dp e 



-pz 



Va{p) 

-Ua{p), 



iFa 2'Kt 

cos 

2 a 



dp e 



\-Ub{p) V + ^l'^ 



where (VaCp); Ua{jp)) and (Vb(p), Ub{p)) are, respectively, the Borel transforms 
of {VAiz'), Ua{z')) and {Vb{z') - z' , Ub{z') - z' - 1/2). Thus the following term 
is added in Ret > 2T 



ir 



a 



m / dpe 



■pz' 



, ( Va{p)\^ ira 2Kt 

' - ' H cos 

2 cj 



--Fe 



-2niz' 



2K^Biz' 



m 



Ua{p)) 
1 



dp e' 



' ( ^b{p) \ 
\-Ub{p)) 



ia 2'Kt 



Sk^Biz 



12 



+ 2K^iB2z' 



2K^B-iZ 



-2-Kiz' 



2K^Bi{t - t2) 



2'Kt. 



l-l) + 2^^*"*2rcos— Stt^Si L 



1 — 

which leads to terms of order of -^e in the outer solution. Therefore this 

1 -azd 1 ^2 

term is negligible when -^e v^<^ <C — e v^<' . 

As a result, the overall contribution of the singular points is the simple sum of the 
contributions from ti, t2, t\, t^- Because ()6-2() has no singularity in the domain 
Re t > 2T, one finally has the following expression of the unstable manifold valid in 
the whole time domain: 



Vu{t)= Vooit) + 0-^1^02 (i) 

+ S{t)Re 



— — e <^ X2{t)e 
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4^1(2) 2nU2 \ (Tik-l? - (1 + fc) 



4/c(fc- 1) 



a;i(t) + X2{t) ) e 



-4yl(l) 



+ 5(t)Re 
+ 5+(t)Re 



-e <^ I y2(i) + (^^TT- I e 



4/1(2) 2..t2 r /rffc- 1)2 - (l + fc) 
e 



4fc(A; - 1) 



yi(t)+a 



_ 2-nit 

e 



(6-5) 



where 5'+(t) = 5(t — 2T) and ^(t) stands for the step function, and A^^\ A^"^^ are 
given previously. Before closing this section, we give an approximate expression 
of the stable manifold Vs{t),Us{t), which is the time-reversal image of the unstable 
manifold: 



Vs{t)= VQQ{t) + CF'^VQ2{t) 

- S{-t)Re 



— : — e " X2{t)e 



5+(-t)Re 



4yl(2) 2^ r/T(fc- 1)2- (1 + A:) 



-e 



4k{k - 1) 



xi{t) + X2it) > e 



us{t)= noo(t) + + a^uo2{t) + a' W - ^^i w) 



5(-i)Re 
5+(-t)Re 



e - y2{t) + a 



_ 2TTit 

e ■y 



4^(2) 



-e 



r(fc- 1)2 - (1 + fc) 
4fc(fc- 1) 



+ y2{t) + o 



yi{t)+a 



y'lit) 



2-Kit 

e " 



(6-6) 



§7. Comparison v^ith Simulations and Reconnection of Unstable 

Manifold 

7.1. Comparison with numerical calculations 

Here we compare the analytical solution H6-5() obtained in the previous section 
with the numerical calculation. (|6-5jl indicates that, when t exceeds t = 0, expo- 
nentially growing oscillatory term appears and, when t exceeds t = 2T, another 
exponentially growing term is added. These terms are exponentially small with re- 
spect to a but they grow exponentially and become dominant for sufficiently large 
t. Thus, added terms describe heteroclinic tangles and a large-amplitude oscillation 




Fig. 4. The analytically constructed unstable manifold (solid line) and numerically portrated phase 
space, (cr = 0.35, k — 0.85) 
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Fig. 5. The analytically constructed unstable manifold (solid line) near (— 7r,0) and time evolution 
of the ensemble (plus-mark), {a = 0.3, k = 0.85) 



appears near (— vr, 0) can be explained with them. This is indeed the case. Fig. 0] 
shows the analytically obtained approximate solution of the unstable manifold from 
(7r,0) (solid curve) and numerically portrayed phase space, and Fig. [Slshows the an- 
alytically obtained approximate solution of the unstable manifold from (tt, 0) (solid 
curve) and the numerically calculated time evolution of an ensemble of 500 points 
on solid curve (plus-mark) near (— 7r,0). As shown in Fig. [51 the analytically ob- 
tained approximate solution (|tj-5|) (solid curve) well reproduces the numerical results 
(plus-mark) , even the large-amplitude oscillation and stretching of the ensemble near 
(— 7r,0). When only the leading order terms of the inner equation are retained, the 
agreement between the analytical and numerical results is not so good (cf . Fig. |B} . 
We also observe that terms of order a extends time range in which the solution well 
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(x10"-3) (x10"-5) 




V-Tl V-Tl 



Fig. 6. The analytically constructed unstable manifold (solid line) near (— 7r,0) and time evolution 
of the ensemble (plus-mark). Only the leading term is taken into consideration, (cr — 0.3, fc — 
0.85) 
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Fig. 7. The analytically constructed unstable manifold. The left figure shows the case of (t = 
0.3, k = 0.3 and the right figure shows the case of o" = 0.3, k = 0.85. 

represents the unstable manifold and, thus, the terms of order a of the inner solution 
are important. 

Fig-El Fig.|Slshow the unstable manifolds when a = 0.30 and k = 0.3, 0.85, 1.0— 
5.0 X 10~^^, respectively. The overall structure of the unstable manifold is similar to 
that of the separatrix except near the fixed point (— vr, 0). Near (— vr, 0), the unstable 
manifold acquires an oscillatory portion and the slope of this portion becomes steeper 
as k increases 1. The behavior of the slope can be understood from the asymptotic 
ratio between the additional terms: 

lim — t\^-^ — = —Vk a . 

t^+oo Vl{t) + TT 2 

In Fig. [71 Fig. jSl one can see a small stochastic domain and this ratio represent an 
unstable direction in the domain near (— vr, 0). 
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Fig. 9. The analytically constructed unstable manifold near (0, vr). The left figure shows the case of 
a = 0.3, k = 1.0-5.0x10"^^ and the right figure shows the case of cr = 0.3, k = 1.0-1.4x10"^^. 



7.2. Reconnection transition 

We remind that the separatrix of the continuous-time hmit equation is the un- 
perturbed unstable/stable manifold. As mentioned in introduction (Fig. the sep- 
aratrix changes its topology depending on the parameter k (reconnection). When 
< < 1, there appears a separatrix connecting (vr, 0) and (— 7r,0). As A; — > 1, the 
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time when the orbit crosses the li-axis diverges and it approaches the separatrix for 
k = 1 connecting (vr,0) with (0, vr). The corresponding topological change occurs for 
unstable/stable manifolds, which will be discussed in subthis section. As shown in 
the previous section, since the approximate stable and unstable manifolds are related 
with each other by a simple symmetry, it is enough to discuss the topological change 
of the unstable manifold. It is interesting to see that, when k is very close to unity 
(i.e., A; = 1.0 - 5.0 x 10""), there appears an additional oscillatory portion in the 
unstable manifold near (0, vr) (cf. Fig. |HJ Fig. EJ. As mentioned before, at /c = 1, 
two segments from (vr, 0) to (0, vr) and from (0, vr) to (— vr, 0) are separatrices of the 
continuous-limit equation. Thus, the perturbed unstable manifold at A; = 1 start- 
ing from (vr, 0) should have an oscillatory portion near (0, vr). The oscillation near 
(0,7r) of the unstable manifold shown in Fig. can be considered as the precursor 
of this oscillation in the unstable manifold at A; = 1. Origins of these behaviors are 
summarized as follows. 

When k is not very close to 1, the unstable manifold reaches near (— vr, 0) before 
additional terms become large. On the other hand, when k is very close to 1, the orbit 
along the unstable manifold takes very long time to reach u-axis because VQo{T(k)) = 
and lim T(k) = oo. Moreover, since the solutions at t = T(k) behave: 



lim , , _ , 

X2{T)\ _ f oo 



lim , . _ , 

fc^i-o \y2[T) J \-oo 

when A; ~ 1 the additinal term is not so small compared to the perturbative solution 

near (0, vr). Therefore the unstable manifold shows a small oscillation near (0,vr). 

We note that the oscillation near (0, vr) is driven by the first additinal term, on 

the contrary, the oscillation near (— 7r,0) is driven by the two additinal terms. We 

further remark that two oscillatory portions in the unstable manifold for A; ~ 1 come 

from the singularities at t = ti , and the singularities t = t2 j ^2 goes to infinity as 

lim Re t2 = oo. In other words, the two sequences of singularities are necessary 
fc— >i— 

for the description of the reconnection of unstable/stable manifolds. 

In Fig. El the unstable manifold for A: = 1.0 — 5.0 x 10"^^ and A; = 1.0 — 1.4 x 
10~ near (0, vr) are shown. Because of the exponential divergence (~ e^*) of the 
additional terms for t — > oo, the approximate unstable manifold H6-5() will not work 
if |ee^*| > 1, i.e. t > f^. When k is very close to 1, T ~ -|ln(l - A;) > 1, and, 
thus, (|H-5jl is not valid for t > T. In this case, (|6-5)) is applicable for t <^T. Indeed, 
when A; = 1.0 - 1.4 x 10"^^ (lESI) shows self crossing curves. Hence, Fig. El is drawn 
by restriction oi t < T. 

§8. Conclusions 



With the aid of the ABAO (the asymptotics beyond all orders) method, we 
have derived analytical approximation of the unstable/stable manifolds, which agree 
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rather well with those obtained by the numerical iteration. The perturbed unstable 
manifold starting from (vr, 0) acquires additional terms which are exponentially small 
with respect to perturbation parameter a, but it exponentially grows with respect 
to time t. Therefore separatrix splitting can be explained by the change of dominant 
terms among perturbative solution and added terms derived by the ABAO method. 
And overall approximation shows a highly oscillatory behavior only near the fixed 
point (— vr, 0) when k is not close to 1. From the comparison with the numerical 
calculation, we check that the approximation breaks dowm when time goes to large 
and that the time regime where the approximation is valid extends thanks to the 
contribution from the higher order terms with respect to a. We also mention that 
when k is very close to 1 we can observe an oscillation near (0, vr), which is considered 
to be the precursor of the heteroclinic tangle in the unstable manifold at A: = 1. In 
this way, even when the heteroclinic tangle exists, the unstable manifold smoothly 
changes its topology as the change of the parameter k. 

In the systems studied so far such as, the standard map, the Henon map, and 
the cubic map, the perturbative solution has only one sequence of singularities in 
the complex time domain. On the contrary, the Harper map exhibits two sequences 
of singularities when k ^ 1. The interference of the contribution from them might 
be expected. We concretely write down the interference term and show that it can 
be negligible. Also we show that it is important, in general, to choose suitable initial 
condition because incorrect choice of the initial conditions lead to Stokes multipli- 
ers(SM) which diverge as A; — > 1. This also shows that even the perturbative solution 
does not have a symmetry with respect to u axis. In a higher order approximation, 
the interference should contribute and such a situation will be discussed elsewhere. 

Before closing this section, we summarize the matching conditions and order of 
the added terms in the ABAO analysis for the following maps (tabled). 

the standard majPE) : f ' ""A = ( ^ 

\Un+l -UnJ \cr Sm Vn J 

the Henon majP-EHJ : f ~ ""A = ( , ,1 

the cubic maiP : ~ ""A = ^ , ,1 

\Un+l-UnJ \-4:aVn{Vn - a){Vn + a) J 

as well as the Harper map. 
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imaginary part of 
a nearby singularity 


matching between 
leading order terms 


order of 
added terms 


Standard map 


J. iTT 

I'C — 2 


vio{tp + (7z) z^e 


1 


Henon map 


J. ZTT 

~ a/2 


<7 ■vio(ip + cz) w z e 


1 V^ir-^ 


Cubic map 


J. iiT 
~ 4a 


n 2-Kltc 

crvioitp + az) ^ z e <y 


1 

^ ^2 


Harper map 


J. ITT 

~ 2Vk 


2-jritc 

^^io(*p + c^z) ze 





Table I. Order of added terms and matching conditions for several maps. The parameter cr plays 
a role of the small parameter, tp stands for the singularity in the upper half plane nearest to 
the real axis and tc — tp — Ketp and vw (t) is the additional terms to the outer equation which 
appears when t exceeds Re tp. 
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Appendix A 

Relation between cases of k > 1 and < k < 1 

In this section, we consider the case of A; > 1. Put 

T = kt, a = ka, 5(t) = v (^—— + a^ , v{t) = u + (^'1) 

then V, u are the solution of the Harper map with parameter \/k: 

v{t + a) — v{t) = — asinn(r) 

u{t + a) — u{t) = — sin'u(r + a) (A-2) 

fv 

Appendix B 

Contribution from complex conjugate singularity 

In the main part of this article, we present details of analytic continuation via 
the behavior of the solution near ti. The purpose of this appendix is to prove that 
the additional term arrising from is just the complex conjugate of that from ti. 
We present the proof for dominant terms, but we only use a real analyticity of 
perturbative solution. Therefore, the same discussion can be used for less singular 
terms. From a real analyticity of {von,uon) and the expansion H2-4|) in a main part, 
{von,uon) has a following expansions near ti and t^ in the complex time domain. 
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Hence, this expansion gives the foUowing asymptotic expansion: 



oo 

y- 



1 = 1 



+ 



oo ^ 



i=l 



(B-l) 



where z and z are defined hy z = {t — ti)/a and z = {t — tl)/a. Therefore, it is 
sufficient to show that the second term generates just the complex conjugate of the 
first term: 



Vadd{z) 

Uadd{z) 



p 



'■0 

(i) 



Let 7' denote symmetric pass of 7 with respect to real axis in p plane. By analitically 
continuing the second term of ()B-1|) from Kez < 0, Imz > to Rez > 0, Imz > 
following term is added. 



dpe-P' Yl 



P 



dpe-P'^ ^' 




i=i ^ ' 



fVaddi} 



\Uadd{z*), 

near t\. Let us suppose that Xadd{t,^) is matched with (|B-2() near ti. Ofcourse, the 
second term is matched with Xadd(t^(^) and it is the complex conjugate of the first 
term, Xaddi't^'^)^ for real t and this fact complete the proof. 



Appendix C 

Borel-trans formed linearized inner equation 



The aim of this appendix is to examine an asymptotic behavior of {Va{z), Ua{z)) 
and {Vb{z), Ub{z)). The Borel transforms {Va{p),Ua{p)) of (Va(z), Ua{z)) and 
{VBip), Usip)) of (VBiz) - z, Ub{z) -{z + \)) satisfy 



where 



{e-'P -l)VA{p) = UA{p)*g{p) 
(1 - eP)UA{p) = Va{p) * f{p) 

{e-P - 1)Vb{p) = Ub{p) * g{p) + ^ff(p) +9'{p) 

(1 - eP)UB{p) = Vb{p) * f{p) + f'{p) 



B 



=«Vbo(z) 



f{p), B 



-iUooiz) 



g{p) 



(C-l) 
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By substituting the power series expansions: 

oo oo 

n=0 n=0 
oo oo -. 

n=0 n=0 

into ()C-1|) and comparing term by term, the coefficients a^, bn are determined re- 
cursively and one has 

AL-i = 0, (Vn) 
AL(-l)'^(27r)2" ^ -i?4n (n ^ oo) 
<-i(-ir+'(2vr)2"-'--i?5 (n-oo) 

A^„(-l)'^(27r)2" ^ -Bin (n ^ oo) (C-2) 

BL-1 = 0, (Vn) 
S^(-l)''(27r)2" ^ Si(2n + 2)(2n + 1) - San (n ^ oo) 
<-i(-ir+'(2^)'""'--^3n (n^oo) 

S2^„(-l)"(27r)2" ^ -5i(2n + 2)(2n + l) -San (n ^ oo) (C-3) 

where Bi = 0.01480, B2 = 0.14, B3 = 0.186, S4 = 3.503, S5 = 5.551. 
Therefore, the following relations hold: 



[Ua{p)) [27r^ B,fi''\p)^ 



-2vrze^--ReSp=2.i (^^j^l'T.) = 2-'^^4. 



-2..e--ReSp... (ll^p;) = 8.^^..^ (_\) + 2.3,^,. Q _ 2.^B,. (j) 

Appendix D 

Borel-transformed first order solution of inner equation 



The aim of this appendix is to examine an asymptotic behavior of the solutions 

4 



Vqi{z), Uqi{z). The Borel transforms Vi{p),Ui{jp) of Vqi{z),Uqi{z) — i^^-r- satisfy 



~ k — 1 k — 1 ~ 

{e-P - l)Vi{p) = i^-g'{p) + i^^g{p) + Ui{p) * g{p) 

(1 - ePpiip) = i^f'{p) + Vr{p) * f{p) (D-l) 
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The power series expansions of Vi{p), Ui{p) is defined by 

Vi{p) = + E , Uiip) = ^^^^ + f: d^p- (D.2) 



n=l 



24 



n=l 



D-l|) and (|D-2|) gives the following form: 



Viip) = ktMp) + ik- i)Vyip) + v,ip) , Ulip) = ktiU.ip) + {k- i)Uyip) + uAp) 

(D-3) 

where Va{p), Ua{p) (a = x,y,z) are independent of k. By substituting ()D-3p into 
(|D4|) . we get 

(e-P-l)y„(p) = ^,(p)*r7 
(1 - enUo^ip) = Va{p) * f 
{a = X, z) 

(e-f - l)Vy{p) = Uy{p) * g + ^g{p) + ^g'{p) 



{l-enUy{p) = Vy{p)*f+lnp) 



By comparing ()D-4|) with HC-1|) . we get 

kti + 1 



Vi{p) = i 

Ui{p)=i 
This estimation gives ()4-33|) . 



24 
kti + 1 
24 



^a(p) + ^Vb{p) 
Ua{p) + ^^Ub{p) 



(D-4) 



(D-5) 



Appendix E 

Choice of initial time and Stokes Multiplier ■ 



In this appendix, the relation between initial time and Stokes multiplier is dis- 
cussed. Here we restrict our discussion to the case of (wo2(i), ^02), which is used to 
obtain our final formula. We start from pertarbative solution of (^02(^)1 1402(0 )• 



1 

V02{t) = 



( sinh ^/kt — ^Jk cosh \fkt ) ( cosh \fkt — \fk sinh \fkt 
c;(t) <^ kt-2Vk^ ' ^ 



(1 + k) cosh^ \fkt — 2\fk sinh \fkt cosh \fkt 



A(k) , , B(k) , , 
+ -r7^a;i(t) + -^X2{t) 



24 



24 



"02 (i) = ^ 



( sinh \fkt — \fk cosh \fkt ) ( cosh \fkt — \fk sinh \fkt 
2y[{t) - yi{t) { kt - 2Vk^ 



(1 + k) cosh^ y/kt — 2y/k sinh y/kt cosh -v/fct 



A(k) , , , , 



24 



24 



34 



S. Ajisaka and S. Tasaki 



We note that {A, B) is chosen as {A, B) = (0, 0) in the main part. In this appendix, 
we put i? = for fo2(— oo) = 0, tio2(— oo) = and leave A. Then, the set of second 
divergent terms in perturabative soluion near ti and t2 are respectively 

~ .( iVk7r-2A{k) + 2 ~ . .^k-l~ \ 
ViiP) = * I ^ yA{p) + -^Vb{p) I 

= M ^ — ^^^P^ + -^^b{p)\ (E-1) 

and 

They derive the following Stokes multipliers. 



Af = -{k - l)7r3B2 - 2A ' 



.R X s 4.kT + i\fk-K - 2Aik) - 2 . , , 

ylf = -(A: - \)tC-B2 + — vr354 (E-3) 

where the upper multiplier is for t\ and the lower one is for t^- We note that the choice 
A()z) = kT = y/k\n^^= makes a symmetric perturbative solution with respect to 
u axis but it also makes the upper multiplier divergent. Therefore this choice should 
be rejected. 

Appendix F 

Matching of higher order Terms 



In this subsection, Matching between inner and outer equation for higher orders 
is discussed. The analytic continuation of a inner solution Vqz to the next sector is: 

Voi{z)+e-^^''Vu{z) 

where Vu, Uu should be determined by the linearized equation. To deal with more 
general case, we analyze the following linearized inner equation: 

oo 

AVi{z,a) = Y,{fin{z)a''V^{z,a) + f2n{z)a^Ui{z,a)} 

n=0 

oo 

AUiiz,a) = Y,{9in{z)a"V^{z,a)+g2niz)a''Ui{z,a)} (F-1) 

n=0 
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Next, we seek to find the solution of the type 



Because ()F-1|) is linear and second difference equation, each ( Yl[^\ 1 has two unde- 



fined constant corresponding to homogeneous solution. It indicates that Vi, Ui has 
the form: 



Ui{z)) [Uu{z) 

oo 



a i=0 



where Aa{o') can be expanded as a power series of a and where 
are defined by the following rules. 

is the indenendent solutions of (Y}^^,^' 



Va{z)\ (Vn{z) 



1. rr / \ is the independent solutions of , ; which satisfy 
\Uaiz)J \Uioiz)J 

lim z"^^Va{z) = l{a = A,B), Ja < Jb 

\z\—>oo 

Note that this condition is not imposed on Ua- 

2. lY}'&^^ is the solution when (Y}'^] = (YA'l) and = 



U?iiz)J \Uw{z)J \Ua{z)J \Uu{z)J \U^,{z) 

{i > / — 1). More concretely, they are inductively defined by the followin; 
equations. 

\U^(^)J fr", \9ijiz)vr_,{z) + g2,{z)W.,{z) 

3. {coefficient of z^-^ and z^b V-^{z)} =0 
In (|2-2|) . boundary conditions for Vni , Uni {n > 1) and e are unsettled. On the other 
hands, Af, Af in (|F-3p is determined by analytic continuation of asymptotic expan- 
sions. Therefore Af, Af determine (|2-2() by matching the asymptotic expansions 
with it. Similar to the case of {V{^, C/f ), solutions of vii{t), vuit) are 

vio\ _ A f VA(t)\ B ( VB(t) 



nioj Wit) J +^0 \uB{t) 



viiit) 



I 

E 

j=0 



A 



1 i-j^ 
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where ( ) denote the homogeneous solution which satisfy: 



1 = hm(t - hy^'^Vait) 



and where Vjj(t), u'^i{t), {I > 1) are defined by following 2 conditions. 
1. f'']^) is the solution when f ""^^^ = f and = (""l^) {0 < i < I - 1) 



2. {coefficient of {t - ti)^" in vfj =0 
The matching condition is: 



e-^^''J2'''^^^^^^ = mY.^'vu{t) (F-4) 



(=0 



As a result, we get 



E{t) = a-^^e-"^^'' 



and the second equality leads 

cf = 0, (j < Jb - Ja) 

4 = '^f-JA-0B^ ij>jB-jA) 

The rest of this appendix is devoted to present a concrete calculation for the Harper 
map. The Harper map has the following values. 

OA = -1, JB = 1 
It indicates the following relation. 

E{t) = = -e^ 

a a 



oo 

i 



Vli{t) 

Unit) 



As noted above, the coefficient Aa{(j) = ""^^anA^ determine whole solutions in 

n=0 

H2-2p . Let us concretely write down the first two terms: 
dpe-P' {Vo{p) + aViip)} 
dp e-P' {Uoip) + aUiip)} 



n 
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A 



B 



{V,^(z) + Vf,{z)a + ---) 
(US,{z) + UR{z)a + ...) 



{y,i{z) + V,i{z)a + ---) 



~2niz 



A^z + a{A^z + ^A^z^) 
A^z + a{Afz-!^A^{z^ + z)) 



-2-Kiz 



where the contour integral is evaluated with the aid of (|4-27)) and (|4-38j) . This leads 
to the following relations 



Af 



-{k 



(F-5) 



The above relations imply two linear relations among Ai, Bi and -B3, which are 
satisfied by the present numerical estimations rather well: 



Ai 



0.99998 ~ 1 



2Ai 
3Bi 



0.99997 ~ 1 



References 

Saito S, Nomura Y, Hirose K and Ichikawa Y H 1997 Separatrix reconnection and periodic 
orbit annihilation in the Harper map Chaos 7 245-253 and references therein. 
Shinohara S 2002 The threshold for global diffusion in the kicked Harper map Phys. Lett. 
A 298 330-334 and references therein. See also S. Shinohara 1999 PhD thesis, Waseda Uni- 
versity. 

3) Ecalle J 1981 Les fonctions resurgentes, 2 (Publ. Math. D'Orsay, Paris) 

4) Lazutkin V F 1984 "Splitting of separatrices for the Chirikov's standard map" VINITI 
6372/84 preprint [Russian] 

Kruskal M D and Segur H 1991 Asymptotics beyond all orders in a model of crystal growth 
Stud.Appl.Math 85 129-181 

6) Lazutkin V F, Schachmannski I G and Tabanov M B 1989 Splitting of separatrices for 
standard and semistandard mappings Physica D 40 235-248 

7) Gelfreich V G, Lazutkin V F and Tabanov M B 1991 Exponentially small splittings in 
Hamiltonian systems Chaos, 1,137- 142 and references therein. 

Lazutkin V F 1991 Exponentially splitting of separatrices and an analytical integral for the 
semistandard map. Preprint, Universite Paris VII 

9) Lazutkin V F 1992 An analytic integral along the separatrix of the semistandard map: 
Existence and an exponential estimate for the distance between the stable and unstable 
separatrices. Algebra i Anahz 4(4) 110-142 (Russian edition). St. Petersburg Math. J. 4(4) 
721-748 (1993) (English edition) 

10) Lazutkin V F, Gelfreich V, and Svanidze N.V 1994 A refined formula for the separatrix 
splitting for the standard map. Physica D 71(2) 82-101 

11) Gelfreich V 1996 Conjugation to a shift and the splitting of invariant manifolds Applica- 
tiones Mathematicae 24,2 127-140 

12) Gelfreich V 1999 A proof of Exponentially Small Transversality of the Separatrices for the 
Standard Map Comm. Math. Phys. 201(1) 155-216 

13) Gelfreich V and Lazutkin V 2001 "Splitting of Separatrices: perturbation theory and 
exponential smallness" Russ. Math. Surv. 56 499-558 [Russian] 

14) Gelfreich V and Sauzin D 2001 Borel summation and splitting of separatrices for the Henon 
map Annal. Instit. Fouriei',51, 513-567 and references therein. 

15) Gelfreich V 2002 Near strongly resonant periodic orbits in a Hamiltonian system Proc. 
Nat. Acad. Sci. USA 99 13975-13979 

16) Gelfreich V and Lerman L 2003 Long-periodic orbits and invariant tori in a singularly 
perturbed Hamiltonian system Physica D 176 125-146 and references therein. 



S. Ajisaka and S. Tasaki 



17) Hakim V and Mallick K 1993 Exponentially small splitting of separatrices, matching in 
the complex plane and Borel summation Nonlinearity 6 57-70 

18) Tovbis A, Tsuchiya M and Jaffe C 1998 Exponential asymptotic expansions and approx- 
imations of the unstable and stable manifolds of singularly perturbed systems with the 
Henon map as an example Chaos 8 665-681; "Exponential asymptotic expansions and ap- 
proximations of the unstable and stable manifolds of the Henon map", preprint (1994). 

19) Nakamura K and Hamada M 1996 Asymptotic expansion of homoclinic structures in a 
symplectic mapping J. Phys. A 29 7315-7327 

20) Nakamura K 1997 Quantum versus Chaos: Questions Emerging from Mesoscopic Cosmos 
(Kluwer, Boston) 

21) Nakamura K and Kushibe H 2000 Beyond- All-Orders Asymptotics and Heteroclinic Struc- 
tures Prog.Theor.Phys. Supplement 139 178-190 

22) Hirata Y, Nozaki K and Konishi T 1999 Exponentially Small Oscillation of 2-Dimensional 
Stable and Unstable Manifolds in 4-Dimensional Symplectic Mappings Prog. Theor. Phys. 

101 No.5 1181-1185 

23) Hirata Y, Nozaki K and Konishi T 1999 The Intersection Angles between N-Dimensional 
Stable and Unstable Manifolds in 2N-Dimensional Symplectic Mappings Prog. Theor. Phys. 

102 No.3 701-706 



